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The time-dependent probability density function of the order parameter of a system evolving 
towards a stationary state exhibits an oscillatory behavior if the eigenvalues of the corresponding 
evolution operator are complex. The frequencies ui n , with which the system reaches its stationary 
state, correspond to the imaginary part of such eigenvalues. If the system (at the stationary state) 
is further driven by a small and oscillating perturbation with a given frequency u, we formally 
prove that the linear response to the probability density function is enhanced when ui = ui n for 
n £ N. We prove that the occurrence of this phenomenon is characteristic of systems that are in 
a non-equilibrium stationary state. In particular we obtain an explicit formula for the frequency- 
dependent mobility in terms of the relaxation to the stationary state of the (unperturbed) probability 
current. We test all these predictions by means of numerical simulations considering an ensemble 
of non-interacting overdamped particles on a tilted periodic potential. 



PACS numbers: 05.40.-a,05.10.Gg,05.70.Ln 



I. INTRODUCTION 

The linear response theory has long been investigated 
for systems at thermodynamic equilibrium Re- 
cently, much effort has been made to develop a parallel 
theory of linear response for systems far from equilib- 
rium and in particular, for those which attain a non- 
equilibrium stationary state (NESS) @. Moreover, the 
study of relaxation phenomena towards a NESS has also 
become an important issue in a variety of systems. Exam- 
ples of the latter are biochemical reaction systems @, [|| , 
growth process [9j and systems in which there occur tran- 
sitions between stationary states (Toj . among others. 

In this work we are concerned with the response to 
a time-dependent perturbation of systems with (contin- 
uous or discrete) markovian dynamics. We report a 
resonant behavior which is exhibited by systems in a 
NESS, and a relationship between such a response func- 
tion and the relaxation properties of the system towards 
the NESS. The occurrence of the resonances in the re- 
sponse function can be intuitively explained as follows. 
The probability density function (PDF) relaxes to the 
stationary state through a superposition of oscillations 
with frequencies {w„}. The latter are given by the imag- 
inary part of the (complex) eigenvalues of a given gener- 
ator C. If this system is further perturbed with an ex- 
ternal forcing oscillating with a given frequency u>, then 
we would expect to have an enhancement of the response 
to the PDF just at oj = oj n . Then, this effect can be 
perceived in some observable quantities like, for exam- 
ple the probability current. That the system must attain 
a NESS is a necessary (but not a sufficient) condition 
for the resonant response to occur. Indeed, it is known 
that for discrete markovian systems the detailed balance 
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(DB) condition implies that the gen erator C can be trans- 
formed into a symmetric one C* [llj . This is then equiva- 
lent to say that DB imply real eigenvalues. Nevertheless, 
for continuous systems it has been derived a symmetry 
property [l2j from the detailed balance, which we show 
leads to a similar result. This excludes the possibility of 
having resonant response for markovian systems fulfilling 
detailed balance. 

At this point it is worth to speak about some other 
works related to the ours. First of all we should point 
out that the resonance phenomenon in the susceptibility 
function was first predicted by Ruelle in [lij [l4| for uni- 
formly and non-uniformly hyperbolic dynamical systems. 
The former systems are known to have the property that 
time averages are uniquely determined by a Sinai-Ruelle- 
Bowen measure. The latter kind of measures have been 
proposed as candidates for non-equilibrium stationary 
states [TEj and some important consequences have been 
inferred from such hypothesis. An example of the lat- 
ter is the Gallavotti-Cohen fluctuation theorem [pel Il7|. 
In the context of dynamical systems, Ruelle stated that 
the susceptibility function should have two types of sin- 
gularities in the w-complex plane: some of them related 
to the unstable dynamics and others corresponding to 
the stable one [13|, [l4[ of the hyperbolic dynamics. Such 
singularities are called resonances (see Ref. [l5[ section 
4.9 and [lj])- In particular Ruelle stated that the sin- 
gularities corresponding to the stable dynamics "would 
correspond to resonances in the 'oscillations of the system 
around its attractor' " (l4) . while the singularities corre- 
sponding to the unstable one are the same as those of 
the spectral density (i.e., the fluctuation-dissipation re- 
lation holds only for the unstable contribution [r| [l4| ) . 
The general framework of Ruelle is applicable to chaotic 
systems and the presence of such resonances has been 
verified analytically and numerically [1 8l — |2Q| ] in specific 
chaotic models. However, in the context of stochastic 
systems attaining a non-equilibrium stationary state such 
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phenomena has not been previously studied (up to the 
knowledge of the author). In this paper we extend part 
of the Ruelle's ideas to the stochastic domain. It is 
worth mentioning that the linear response for the case 
of stochastic systems is very different from that of the 
chaotic case [fjj. For example, difficulties related with 
the possibility that the Sinai-Ruelle-Bowen measures be 
singular (i.e. not absolutely continuous with respect to 
Lesbegue) and the decomposition of the state space into 
stable and unstable directions, are absent in our case. 
Thus, the phenomena observed in the linear response for 
the models analyzed in [181 - 120} are different and more 
complex than the behavior found in the systems explored 
in this work. 

This paper is organized as follows, in Section |n] we 
prove that the resonance in the linear response may oc- 
cur if the system reaches a NESS via oscillations due 
to complex eigenvalues. In Section IIIII we show the oc- 
currence of the resonant phenomenon in titled periodic 
potentials by means of numerical experiments. Finally 
in Section lIVI we give our conclusions. Some appendices 
are included containing detailed calculations. 



II. THE LINEAR RESPONSE AND THE 
RELAXATION OF THE PDF 



now follow the usual way to obtain the linear response to 
the PDF [6| . First we assume that the perturbed solution 
can be written as 
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In the above ansatz, the deviation from the stationary 
state of the perturbed system, with respect to the un- 
perturbed one, is represented by e\R)e lult . As \p p ) is a 
normalized probability vector, it follows that \R) is a vec- 
tor parallel to the simplex Ax, i.e., (— \R) — 0. We will 
refer to \R) as the response to the stationary PDF. It 
is clear that in the ansatz Q we use the basic assump- 
tion that the perturbed PDF will eventually be periodic, 
oscillating with the same frequency of the perturbation. 
Using Eq. ^ in §5$ we obtain the following equation for 



(iu-£)\R) = C p \P ). 



(5) 



The linear response to the PDF has the formal solution 
\R) = (iuj — £ p \Po). We now use the fact that the 
inverse operator (s — C) can be formally written in an 
integral representation as [2lj |. 
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Consider an unperturbed markovian system with state 
space X, whose dynamics is defined by a generator C. 
Throughout this work we will use the Dirac's bracket 
notation as in Ref. 12]. We denote the time-dependent 
PDF of such a system by p(x,t) = (x\p(t)), with x 6 X. 
Then, the evolution equation for \p) can be written as, 



d\p) n , 

By using the separation ansatz \p{t)) 
at the eigenvalue problem for C, 

C\P X ) = X\P X ). 



(1) 

e At |P a), we arrive 
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It is known that the largest eigenvalue is A = 0. As- 
suming non-degeneracy of the latter, the corresponding 
eigenvector \P ) is the unique stationary solution of (fTJ) 
on the set Ax of all normalized PDFs. The set Ax 
of probability densities is commonly called the simplex 
space and is defined as Ax := {\p) '■ {—\p) = !}■ Here 
|— ) is a probability density (called the flat distribution) 
perpendicular to the simplex Ax and is defined by the 
equation (— \x) — 1. It is clear that all the other eigen- 
values satisfy Re(A) < 0. Moreover, the eigenvectors cor- 
responding to the non-vanishing eigenvalues are parallel 
to the simplex Ax, i.e., (— \P\) = for all A ^ 0. 
We now consider the perturbed evolution equation, 



d\p P ) 
dt 
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where e > is a small control parameter, w > is a 
given frequency and C p is a perturbing operator. We 



for Re(s) > 0. If we set s = S — iu, with S > 0, we can 
write R as, 
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At this point it is necessary to point out that the perturb- 
ing operator is such that (— \C P — 0. The latter follows 
from the probability conservation of both, the perturbed 
and the unperturbed system [l2j |. Such a property lets 
us to interpret e ct £ p \Po) as the time-evolution of a vec- 
tor parallel to the simplex with initial condition C p \Pq). 
With this in mind we can see that e C p Pq can be ex- 
panded in a series of the eigenvectors of £ as, 
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Notice that the eigenvector corresponding to A = is 
excluded, since e ct C p \Po) has no components perpendic- 
ular to the simplex. If we put the last expression into © 
we obtain, 
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which results in, 
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The mechanism leading to the resonant response becomes 
clear from Eq. ([5]) : the enhancement of the amplitude of 
the response to the PDF, as a function of to, is a con- 
sequence of the oscillations of the PDF in the simplex 
space towards its NESS. At this point it is important to 
stress that under our assumptions, the resonant behavior 
predicted by Eq. |8| is a characteristic of non- equilibrium 
stationary states which is not present in systems at equi- 
librium. This follows from the fact that DB implies that 
the eigenvalues of C are real. The latter is a consequence 
of the fact that DB can be used to transform the operator 
C in a symmetric one. This has been shown in Ref. 
if £ is a transition matrix with a discrete state space X. 
In Appendix |A] we give a proof, based on a symmetry 
property derived from the DB by Kurchan [12j . for the 
equivalent statement when C stands for a Fokker-Planck 
operator. 

The above-described resonant behavior occurs at the 
level of probability densities since it corresponds to an 
enhancement of the PDF oscillations. However, in a real 
experiment it is not always possible to measure the PDF 
itself. Then, it is important to explore how this behavior 
is displayed by a given observable. Here we analyze the 
case in which such an observable is the probability cur- 
rent. Consider the perturbed system ([3]) in the asymp- 
totic (time-dependent) state \p p ) and let JC(t) be the op- 
erator for the probability current. The expected value of 
the operator K{— oo) gives the probability current at the 
stationary state, 

jo = (-|/C(-oo)|P ). 

In a number of cases the generator C can be factorized 
as £. = VK-o- For example, if £ represents the Fokker- 
Planck operator, V — d x and /Co = / — P~ 1 d x , with / as 
minus the gradient of a given potential V. In a similar 
way we also assume that the perturbing operator C p can 
be factorized as T>1C P , and therefore, the probability cur- 
rent operator can be written as K.(t) = K.Q + elC p e tuJt H(t). 
Here H(t) stands for the heavyside function. From the 
form of K(t) we infer that JC(— oo) = /Co- The probabil- 
ity current for the perturbed system j p can be obtained 
by 

Jp (t) = H(K +£/C p e^ff(i))Mt)). 

With the above expression we can calculate the proba- 
bility current at first order in e. Using Eq. (j4|) we obtain 
for t > 0, 

J P (t) = Jo + ee iut ((-\JC P \P ) + {-\Ko\R(u))) , 

and from the last expression we identify 

MM :={-\JC p \P ) + {-\)C \R(lj)), (9) 

with the linear response to the probability current (the 
frequency dependent mobility). Substituting \R) by the 
expression given by Eq. ([6]) into the above equation we 
obtain 

/>oo 

MM = (-\)C P \Po) + lim / dte- st (~\ICoe ct C p \P }, 



or equivalently, 

/>OC 

H(w) = (-\Kp\P ) + / dte^(-|/Coe £i £ p |Po). (10) 
Jo 

In order to give some physical interpretation to the above 
expression lets consider again the evolution equation ([T]) 
of the unperturbed system. First, we prepare the system 
with the initial PDF |p*(0)> := \P ) + "fC p \P ), where 
7 is, in some extent, arbitrarily chosen, just requiring 
that \p*(0)) > 0. Here \Pq) represents the stationary 
state of our system and jC p \P ) is a vector parallel to 
the simplex. This observation makes it clear that |p*(0)) 
fulfill the normalization condition. The formal expression 
for the time-dependent PDF, with the described initial 
condition, is given by \p*(t)) = \Po) + e ct C p \Po). If we 
use the latter to calculate the probability current j*(t) 
for this specific system we found that, 

f(t) =jo+jHlCoe ct £ p \Po). (11) 
The last results lets us rewrite Eq. (fTOl) as, 

f i(u) = (-\IC p \Po)+ 1 - 1 dte iut (f(t)-j ), (12) 
Jo 

which relates the relaxation towards the stationary state 
with the linear response to an external perturbation. 



III. RESONANT RESPONSE IN TILTED 
PERIODIC POTENTIALS: NUMERICAL 
EXPERIMENTS 

In this section we show the presence of the resonant 
behavior predicted by Eq. (|5]) in tilted periodic poten- 
tials by means of numerical simulations. In particular we 
test the relationship between the relaxation to the steady 
state and the response to an external perturbation given 
by Eq. (1121) . First we consider an unperturbed system 
consisting of non-interacting overdamped particles in a 
tilted periodic potential V(x) with period L and tilt Fq 
(i.e. V(x + L) = V(x) — FqL) subjected to a gaussian 
white noise. The stochastic dynamics of these particles 
is ruled by the Langevin equation, 

where 7 is the friction coefficient which, as usual, will be 
set to one. The function f(x) is minus the gradient of 
the tilted periodic potential V(x). With this we have that 
f(x) is periodic with period L. The term represents 
a gaussian white noise with zero mean and correlation 
(£(0)£(t)) = 2/3 _1 <5(t). Here ft is the inverse temperature 
of the system. 

The Fokker-Planck equation corresponding to the 
Langevin dynamics (1131) is given by, 

% = £fpP, (14) 
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where Cpp is given by 



(15) 



Unfortunately the exact eigenvalues A„ for the Fokker- 
Planck operator £pp are not generally known. Never- 
theless, the author and coworkers have given an approx- 
imation to A„ in [2^ | . Such an approximation is valid for 
large values of the tilt, F » max{/(x)} — min{/(x)}, 
and is given by A„ = i2i:nj / L+{2irn/ L) 2 D e Q. Here jo is 
the probability current, D e ff is the effective diffusion co- 
efficient, L is the period of the (tilted) periodic potential 
and n£N. From this it follows that w n = 27m jo / L is an 
approximation to the natural frequencies of the system. 

According to Section UH we have that if this system 
is perturbed by a time dependent periodic forcing sF(t), 
characterized by a frequency w, then, a resonant behavior 
in the PDF will be exhibited at some u> = u>*. Such a 
perturbed system is modeled by the Langevin equation, 



7§ =/(*)+ 



(16) 



with 7 set to one. To achieve the numerical experiments 
showing the occurrence of resonances at the frequencies 
given above, we chose the tilted periodic potential as 
V{x) = cos(27ra)/(27r) — xF . With this potential we 
have that f{x) = sin(27ra;) + Fq and that L = 1. The 
external forcing F(t) is taken as F(t) — cos(wi) and e, 
the small control parameter, is fixed to e = 0.05. 

In Fig. [IJi we show the mobility (the amplitude of the 
current oscillations with respect to the stationary cur- 
rent) as a function of the frequency for the deterministic 
tilted periodic potential with tilt F = 1.80. In such 
a figure we can observe that at zero temperature the 
system exhibit several resonant peaks at the frequencies 
given by the formula ui n = 2-KnvjL (see Eq. (|B7[) from 
Appendix [B]). Here the mean velocity is v « 1.5, and 
therefore uj n « 9.4n. This behavior can be readily ap- 
preciated in Fig. [TJi. In Fig. [TJ) we show the behavior for 
the system at a finite temperature chosen as = 0.01. 
In this figure we see that frequency-dependent mobility 
/i(w) (filled circles) displays a behavior similar to that of 
its deterministic counterpart. It exhibits a resonant peak 
at the natural frequency (which is estimated, according 
to Ref 22] as wo = 2-Kjo/L « 9.4). However, the succes- 
sive resonant peaks seems to be suppressed by the noise 
in the thermalized system. 

Besides the response to the probability current, in 
Fig. [TJd we have also plotted the numerically obtained 
"relaxation curve" (solid line), i.e., the right-hand side 
of Eq. (TT21 . Within the accuracy of our simulations we 
note that a good approximation of the relaxation curve 
to the response to the probability current at the station- 
ary state given by Eq. (TT21) . We now explain how the 
"relaxation curve" was calculated. First we prepared 
the unperturbed system with the initial PDF given by 
-Pinitiai(^) = Pq{x) + j£pPo(%) with <y = 0.08. Here the 
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FIG. 1. (Color online) Frequency-dependent mobility h(uj) 
for the tilted periodic potential V(x) — cos(27rx)/(27r) — xFo 
with Fq — 1.80 (a) at zero temperature and (b) at a tem- 
perature = 0.01. The external perturbation is taken as 
0.05 x cos(ojt) (i.e. e — 0.050). In (a) we can see that the de- 
terministic response exhibits several peaks corresponding to 
the resonance frequencies u n = 2nnv/L (see text). The inset 
shows an amplification in the region near the main resonance 
peak occurring at uj = 9.4 to better appreciate the complex 
structure of the response. In (b) we observe a similar be- 
havior of the response vs. frequency at a finite temperature. 
However, most of the resonant peaks observed in the deter- 
ministic case disappear. We also see that the relaxation of 
the probability current obtained by (|17p (solid line), and the 
response to the stationary current (filled circles) are related 
according to the expression (|12|) (see text for details). 



perturbing operator C p is given by C p — — d x and Po(x) 
is the stationary PDF takes the form, 



Po(x) 



1 

An 



exp[-pV(x)] 



x+L 



exp[/3V(y)]dy, 



with An a constant such that j Q Po(x)dx = 1. Once we 
have prepared the system in such state, we let it evolve 
to the stationary state Po(x). We then measure the time- 
dependent probability current j re iax(i) from t — up to 
a given time T which is large enough to assume that the 
system has reached the stationary state. In Fig. [5^ we 
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FIG. 2. (Color online) The relaxation of the probability cur- 
rent towards the stationary state. Figure (a) shows the prob- 
ability current jreiax as a function of time for the unperturbed 
system (O with V(x) = cos(2ttx)/ (2n) - xF , F = 1.80 and 
/3 _1 = 0.01. We obtained this time series through a simula- 
tion of 20000 particles obeying the Langevin equation (|13[) up 
to a time t — 100 arb. units. The particles were initially dis- 
tributed according to the PDF P (x) + 'yP' (x) with 7 = 0.080. 
(b) The real part of the Fourier transform of j' re iax — jo- To 
calculate the Fourier transform we used a time series for the 
probability current shown in (a). After that we used for- 
mula (|17|) to obtain the "relaxation curve" shown in Fig. [Tb 
(solid line). We proceeded in the same way to calculate the 
"relaxation curves" shown in the subsequent figures. 



plotted jreiax (i) for the potential V(x) and the parameter 
values used to perform the numerical experiment corre- 
sponding to Fig. [T}d. As stated by the right-hand side of 
the equation (fF2|) . we use the time series obtained above 
to calculate the real part of the Fourier Transform of 
Jreiax (t) — jo shown in Fig. (2b- This quantity is then used 
to calculate 

/>oo 

(-|/C p |P )+7 _1 / e^OrelaxW-jo)^ (17) 

Jo 

which corresponds to the "relaxation curve" displayed in 
Fig-IU 3 (solid line). Note that the first term in the above 
expression equals one since JC P = Id in this case. 



FIG. 3. (Color online) Frequency-dependent mobility n(co) 
for the tilted periodic potential V(x) — cos(2ivx)/(2n) — xFo 
with Fo — 0.90 (a) at zero temperature and (b) at a tem- 
perature = 0.01. The external perturbation is taken as 
0.05 x cos(ojt) (i.e. e = 0.050). In (a) we can see that the de- 
terministic response is monotonically decreasing curve. In (b) 
we observe a non-monotonous behavior of the linear response 
(filled dots), which is quite different from its deterministic 
counterpart. We also compare the relaxation of the current 
given by (|17[) (solid line) with the linear response (filled cir- 
cles) in order to test the relaxation-response relation given 
by CCU). 



In Fig. we show the mobility as a function of the 
frequency for the deterministic tilted periodic potential 
for a tilt Fo = 0.90. The behavior of the amplitude os- 
cillations in this case is a monotonically decreasing curve 
whose minimum occurs at uj = 0. This behavior is con- 
sistent with the analytically predicted in the Appendix |B1 
for this case. In Fig. [3Jd we show the mobility for the very 
system but at a temperature = 0.01. All the other 
parameters remain the same as in the case F = 1.80. 
In this figure we can see that the response becomes a 
non-monotonous curve. This is an important feature of 
the tilted periodic potential below the critical tilt since 
the noise induces a behavior which is quite different from 
its deterministic counterpart. This non-monotonous be- 
havior can be seen as emerged from a competition be- 
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FIG. 4. (Color online) Frequency-dependent mobility fi(u) for 
the tilted periodic potential V(x) = cos(27ra;)/(27r) — xFq at 
temperature /3 _1 = 0.01 (a) for Fq — 0.98 (below the critical 
force F c = 1) (b) for Fo = 1.50 (above the critical force F c — 
1). The external perturbation is taken as 0.05 x cos(tJi) (i.e. 
e = 0.050). The graphs compares the mobility obtained by 
the response of the system (filled circles) and the one obtained 
by the relaxation of the probability current (solid line) given 
by (|12|) . The latter formula is shown to be right within the 
accuracy of our numerical simulations. 



tween the deterministic dynamics and the noisy escape 
rate. On the one hand, at low frequencies we have that 
the dynamics of the system is dominated mainly by the 
noise. This is due to the fact that the slow variation of the 
perturbation lets the particles escape from one potential 
well to another in every time period. Thus, many par- 
ticles increase and decrease their velocities, giving large 
amplitudes to the current oscillations. As the frequency 
increases, the number of particles escaping from the po- 
tential wells during each oscillation diminishes. Then, the 
amplitude of these current oscillations no longer profits 
from the escape events and the amplitude of the current 
starts to behave as in the determinist case. This behavior 
can be observed in Fig. [3}d: above some frequency value 
(approximately at u> = 5) the response curve mimics the 
form of its deterministic counterpart shown in Fig. 
Figure [HJs also shows that the relation between the relax- 



FIG. 5. (Color online) The Fourier transform of the time- 
dependent probability current towards the stationary state 
for (a) F = 0.90 and (b) F = 1.50. In both cases 
we prepared the system at three different initial states: (i) 
Finitiai = Po(x) + jPo{x) with 7 = 0.08 (black line), (ii) 
f initial = S(x) (red line) and (iii) P in iti a i = 1/L (green line). 
None of the two latter initial conditions reproduces the be- 
havior of response to the probability current. Nevertheless, 
in all cases the curves obtained show a resonance-like behav- 
ior near the location of the "true resonance" . Thus, we can 
say that these experiments are able to qualitatively indicate 
where the resonance peaks are probably located. 



ation and the response (TTS]) also holds for this case within 
the accuracy of our numerical experiments. 

In Fig. [H we test again the relaxation-response re- 
lation given by (fT2j) for other values of the tilt. In 
Figs. Hi and gfc> we fix F = 0.98 and F = 1.50 re- 
spectively. For such values of the tilt we measured the 
probability current j'o(-Fb) f° r tb unperturbed system giv- 
ing jo (0.98) = 0.28 ± 0.005 and j (1.50) = 1.12 ± 0.005. 
According to Ref. [22| this gives the "natural frequencies" 
w (0.98) rj 1.74 and w (l-50) 7.0 respectively. We per- 
turb the system with the same time-dependent forcing 
used in the experiment described above. We then cal- 
culated the mobility shown in the mentioned figures by 
means of the linear response to the perturbation (filled 
circles) and by means of the relaxation of the probability 
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current (solid line). We observe again a good agreement 
of expression (|12p in both numerical experiments. 

In the case Fq = 0.98 it is interesting to notice that 
although the response is enhanced at the predicted fre- 
quency wo(0.98) ~ 1.74, it does not seem to correspond 
to a maximum located at such frequency. Rather, we 
observe that the maximum seems to occur at lo — 0. Ac- 
tually, in our simulation the response is measured for w's 
bigger than one. However we are not able to observe any 
maximum in the interval 1 < u < 2. This could probably 
be due to either the lack of statistics or the fact that the 
imaginary part of the (non-vanishing) largest eigenvalue 
of such a system is actually smaller than the predicted 
in H|. 

On the other hand, in the numerical experiments per- 
formed, we have seen that Eq. (fl"2"|) is right within the 
accuracy of our simulations. The expression ([12]) states 
that the relaxation from a given initial state to the NESS 
is related to the response to the stationary PDF. The 
initial state referred above is very specific and it is in- 
teresting to see what happens with different initial con- 
ditions. In Fig. [5^ and [5}d we calculated the "relaxation 
curves" for the cases F — 0.90 and F — 1.50 respec- 
tively. In both figure the initial conditions used were (i) 
initial = Po(x) + iPq{x) with 7 = 0.08 (black line), (ii) 
^initial = S(x) (red line) and (hi) P in itiai = l/L (green 
line). It is possible to observe that, although the initial 
conditions (ii) and (iii) do not reproduce the behavior of 
the response, they do indicate qualitatively the position 
of the resonance peaks. This fact is important since it 
can be used to predict the resonance character of a given 
system, not by performing directly the perturbation on 
it but observing its relaxation properties. 



IV. DISCUSSION AND CONCLUSIONS. 

We have shown the existence of a resonant behav- 
ior that is exclusive of non-equilibrium stationary states. 
The requirement that the system be in a NESS is a nec- 
essary but not a sufficient condition for the resonant be- 
havior be present. We proved the latter statement for 
systems with markovian dynamics which includes dis- 
crete and continuous systems. In particular, we have also 
found that the frequency-dependent mobility of a given 
system is related to the relaxation of the probability cur- 
rent to the stationary state with a specific initial con- 
dition. We verified the presence of these resonances for 
the unperturbed system which consists of an ensemble of 
non-interacting overdamped particles in a tilted periodic 
potential. It has been shown that the correspondence es- 
tablished between the mobility and the Fourier transform 
of the relaxation of probability current (Eq. (TT2l ) agreed 
with our numerical experiments for several values of the 
strength of tilt. Finally, we analyzed the relaxation prop- 
erties of the system for other initial conditions different 
from the "specific" mentioned in the text. We found that 
the relaxation curves obtained do not reproduce the re- 



sponse of the system, but they do exhibit a "peak" near 
the resonant frequencies. On this basis, we say that it 
could be possible to predict the resonant behavior of a 
given system by analyzing its relaxation to the NESS. 

We would like to remark that our results could be use- 
ful in understanding some phenomena like the relaxation 
properties of systems that attain a NESS 0413 or the 
so-called resonant transport [23, [24| . It would also be in- 
teresting to see how the relaxation- response relation (Ti"2l , 
is connected to noise-enhanced stability of metastable 
states [H423|. 
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Appendix A: Detailed balance implies real 
eigenvalues 

Let the elements of T be denoted, as usual, by Tij, 
with i,j € X. The detailed balance condition states that 
Tijfti = Tj^nj, where 7r stands for the stationary proba- 
bility vector for T. Following Ref. [Tl| . we should notice 
that the latter can also be written as n~ 1 ^ 2 Ti ^ X J 2 = 

n i 1 ^ 2 ^j',i 7r j^ 2 - ^ we define T* } :— it - 1 ^ 2 T 1 ^tt 1 J 2 we can 
see that the matrix T* is symmetric and therefore its 
eigenvalues are real. Additionally, we can see that T* 
is obtained by a similarity transformation from T, since 
T* = S~ 1 TS where, S% t j = Ttj This implies that 

T* and T have the same eigenvalues. Then it follows that 
if T comply with the detailed balance condition, then it 
has only real eigenvalues. The proof of the latter state- 
ment for C as a Fokker-Planck operator is carried out in 
a similar way. First assume that C comply the detailed 
balance condition and define Pq(x) := (x\Pq) as the cor- 
responding equilibrium stationary state. Then note that 
the transformation defined by S = Pq(x) 1 ^ 2 \x)(x\ has 
inverse S^ 1 = Pq(x)^ 1 > 2 \x}(x\ since, 

SS- 1 = (p (y) 1/2 b)(y|) (P (x)-V 2 \x){x\), 

= Po(y) 1 / 2 Po(x)-^ 2 S(x,y)\y)(x\. 
= Id, 

where Id stands for the identity operator. On the other 
hand we have that 

5 2 = (p (y) 1/2 \y)(y^(Po(x) 1/2 \x)(x\y 

^P (y)^ 2 P (x)^ 2 S(x,y)\ y }(x\. 
= P {x)\x)(x\, 
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In Ref. [12J it has been shown that the detailed balance 
condition leads to the following symmetry property for 
the Fokker-Planck operator, 

QT l CQ = C\ (Al) 

where Q is defined thorough the equation, 

Q\x) =P Q (x)\x). 

It is clear that a representation for Q is given by S 2 . This 
makes it possible rewrite Eq. (|A1[) as, 

(ssy 1 css = tf, 

or, equivalently, as, 



(s- l cs 



t 



(£*) f 



(A2) 



The above equation states that, if the balance condition 
holds, there is a similarity transformation making the 
Fokker-Planck operator hermitian. Then, detailed bal- 
ance condition implies that all the eigenvalues of C are 
real. This proves that a system having a equilibrium 
state does not exhibit the resonant behavior predicted 
by ©. 



Appendix B: Tilted periodic potentials at zero 
temperature 

Consider an ensemble of non-interacting overdampcd 
particles in a tilted periodic potential V(x) — V p (x) — 
Fqx. Here V p (x) represents the periodic part of the po- 
tential, which has period L, and F$ represent the tilt of 
the potential. The deterministic evolution equation for 
one overdamped particle in such a potential is given by, 



(Bl) 



where f(x) is the gradient of the tilted periodic potential 
and f p (x) = —Vp(x) is the gradient of the corresponding 
periodic part. We will first consider the case in which Fq 
exceeds the so called "critical tilt" F c := max{|/ p (x)|}. 
In such a case Eq. (|B1[) have only running solutions and, 
as we anticipated in Section IIIII it exhibits a resonant 
behavior in the frequency-dependent mobility. 

Let x(t) be a solution of (IB1[) with a given initial con- 
dition xq and consider the perturbed system 



dy 
dt 



f(y) + eF{t) 



(B2) 



where e < 1 is a dimensionless control parameter and 
F(t) is a time-dependent periodic external forcing with 
period T. For e small we can measure the response of the 
system as the deviation of the perturbed trajectory y(t) 
from the unperturbed one x(t) having the same initial 
condition y(0) — yo — xq. Then we will analyze the the 



behavior of the "linear response" r/(t) :— (y(t) — x(t))/e. 
From the above we can see that the evolution equation 
for r)(t) is approximately given by, 



dr) 
It 



= f(x(t))r 1 + F(t), 



(B3) 



at first order in e. The above equation can be easily 
solved if we notice that the integration factor v(t) is given 
by, 



:= exp (- J f 



(x(s))ds 



(B4) 



This lets us write a solution for r)(t) as follows, 

1® = / »(s)F(s)ds 
f (*) Jo 

The integration factor v can be written in a more com- 
pact form if we use the fact that the solution x(s) is 
always increasing and therefore invertible. A change of 
variables in the integration in Eq. (|B4I) shows that, 



u{t) 



/(*(*))' 

which can be used to express r\ as, 



V(t) = /(*(*)) 



o /(*(«))' 



ds. 



(B5) 



In Appendix [C]we show that any running solution x(t) 
of the unperturbed system with Fo > F c has the property 
x(t + t) = x(t) + L, where r is the time that the particle 
takes to travel one spatial period. From this property it 
follows that f(x(s)) is periodic with "natural frequency" 
luq = 27t/t, which is proportional to the mean velocity 
v = L/t. If we denote by c„ the coefficients of the Fourier 
series of l/f(x(s)), and take F(s) = e tu}t , we obtain for 
V, 

e i{w—nu a )t _ i 
n{t) = f(x(t)) Cn • (B6) 

n 

From this we can see that the deviation of the perturbed 
trajectories from the unperturbed ones is enhanced (and 
even diverges) when the frequency of the driving force is 
2-7r/L times the mean velocity of the particle, 



2nv 



(B7) 



To tackle the case Fo < F c we should notice that 
Eq. (|B3| and its solution (|B5|) is valid in this case. The 
only difference is that x(t) has no longer the property 
x{t + t) = x{t) + L. Instead, the solution x(t) to the 
differential equation (IB1|) has the asymptotic behavior 
x(t) = x* + Co cxp(— Xt) for t — > oo. Here x* is a root of 
f(x), X = f'(x*) and C = x Q - x* . Thus, for large t we 
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can expand f(x(s)) and f(x{t)) around x* to obtain for 
V, 

r){t) = f e~^ t - s '>F(s)ds. 
Jo 

To obtain the linear response to the particle current we 
fix F(t) — e %lJjy and take the time derivative to r\ giving, 



m = i 



A + iui 



A + ius 



Prom the latter we see that f] is asymptotically periodic, 
and the amplitude of (the real part of) such a function 
is given by, 

A(u) 



A 2 + oj 2 ' 

which is a monotonically increasing function of uj. 

Appendix C: Analytic running solutions. 

The evolution equation for overdamped particles in a 
tilted periodic potential V(x) = V p — xF is 



dx 
~dt 



/(*) 



(CI) 



where f(x) = f p (x) + Fq is the gradient of V(x). Since it 
is assumed that Fq is greater than the critical tilt F c := 
max{|/ p (a;)|}, we have that f(x) has no real roots. This 



means that we can integrate the above equation with the 
initial value x(0) — x Ql giving, 



x(t) 



dy 



o m 



(C2) 



Define t as the time that the particle takes to travel 
one spatial period, i.e., 



dx 



o fix)' 



and notice that this time does not depend on where the 
particle starts its motion. Now define the invertible func- 
tion 



h{x) 



L f x dy 



o f(y)' 



Due to periodicity of f(x) we have that h(x+L) = h(x) + 
L. In terms of this function we can write Eq. (|C2[) as, 



Lt 



h{x(t)) - h(x ). 



The invertibility of h lets us write x(t) as, 



LA 



(t) = h- 1 [ h(x ) + —j 



(C3) 



It is easy to see that h 1 has also the property h 1 (x 
L) = hr 1 {x) + L, from which we obtain that x(t + r) 
x(t) + L. 
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